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LOFAR: Detecting Cosmic Rays with a Radio Telescope 
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Abstract: LOFAR (the Low Frequency Array), a distributed digital radio telescope with stations in the Netherlands, 
Germany, France, Sweden, and the United Kingdom, is designed to enable full-sky monitoring of transient radio sources. 
These capabilities are ideal for the detection of broadband radio pulses generated in cosmic ray air showers. The core 
of LOFAR consists of 24 stations within 4 square kilometers, and each station contains 96 low-band antennas and 48 
high-band antennas. This dense instrumentation will allow detailed studies of the lateral distribution of the radio signal 
in a frequency range of 10-250 MHz. Such studies are key to understanding the various radio emission mechanisms 
within the air shower, as well as for determining the potential of the radio technique for primary particle identification. 
We present the status of the LOFAR cosmic ray program, including the station design and hardware, the triggering and 
filtering schemes, and our initial observations of cosmic-ray-induced radio pulses. 
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1 Introduction 

Radio emission from cosmic ray air showers was discov- 
ered in 1965 by Jelley and collaborators [1]. Advances in 
high-speed digital electronics have expanded the options 
for filtering, triggering, and recording the radio signal and 
have led to renewed interest in the technique. This win- 
dow on high-energy cosmic rays has a duty cycle of nearly 
100% and can in principle, like the air fluorescence tech- 
nique, provide information on shower development in the 
atmosphere. 

The LOFAR Prototype Station (LOPES) successfully de- 
tected the coherent radio pulses from air showers with en- 
ergies from 5 X IQi^ to 5 x IQi"^ eV in a frequency range 
of 43-73 MHz. LOPES verified the primary nature of the 
emission from the air shower: the interaction of shower 
electrons and positrons with the Earth's magnetic field [3]. 
Since then, significant progress has also been made in the 
theoretical understanding of the radio signals [4], with cur- 
rent work focused on sub-dominant emission mechanisms, 
such as charge excess and Cherenkov contributions. 

LOFAR, the Low Frequency Array, is a digital radio tele- 
scope with a dense core of antennas capable of detecting 
air shower pulses and recording the signal polarization and 



lateral distribution with unprecedented accuracy. Under- 
standing and modeling all details of the radio emission will 
allow a full demonstration of the potential for determining 
fundamental cosmic ray properties from the radio signal, 
such as primary particle type [5, 6]. 

2 LOFAR 

LOFAR is a digital aperture-synthesis radio telescope con- 
sisting of at least 48 stations. 24 core stations are concen- 
trated within a 4-km^ area near Exloo in the Netherlands, 
16 remote stations are distributed up to 65 km from the 
core throughout the northern Netherlands, and 8 interna- 
tional stations are spread between England, France, Ger- 
many, and Sweden. Unlike typical ground-based cosmic 
ray detectors, the density of LOFAR stations from the core 
falls off smoothly, from six stations in the "superterp" (the 
360 m diameter island at the center of the LOFAR core; 
see Fig. 1), to effectively isolated remote and international 
stations. 

Each station consists of either one or two high-band an- 
tenna (HBA) fields, covering the frequency range 110-250 
MHz, and a low-band antenna (LEA) field, covering the 
10-80 MHz range. The HBA fields consist of between 24 
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Figure 1 : Photograph of the LOFAR superterp - pictured 
are twelve HBA core fields (grey squares) and six LBA 
fields (barely visible). 



Station Type 


Number 


LBA 


HBA 


Core 


24 


96* 


2x24 


Remote 


16 


96* 


48 


International 


8 


96 


96 



Table 1: Summary of the LOFAR station types. *Only 48 
can be used simultaneously. 

and 96 tiles of 4 x 4 bowtie antennas, while the LBA fields 
each consist of 96 dual-polarization droop-dipole antennas, 
48 of which can be used simultaneously (on the core and 
remote stations). The LOFAR Cosmic Rays Key Science 
Project aims to use the LBAs to detect atmospheric cosmic 
ray air showers in parallel with standard LOFAR opera- 
tions. 

As a fully digital telescope, LOFAR uses digital signal pro- 
cessing to focus the signals from individual station ele- 
ments, in contrast to traditional telescopes which use ge- 
ometrical optics (curved mirrors). Signals from individ- 
ual elements (LBA dipoles or HBA tiles) are digitized at 
200 MHz' with 12-bit precision. The digitized signals 
are divided into 512 frequency sub-bands, up to 240 of 
which can be added together with appropriate phase off- 
sets to form up to eight station beams viewing between 
^ 2 and ^ 1700 square degrees each, depending on the 
frequency selection and antenna configuration [2]. These 
beams are then sent to the LOFAR central processing hub 
in Groningen for cross-correlation (imaging) or further 
beam-forming (e.g. for pulsar searches). 

In parallel to station beam-forming, the digitized time- 
domain signals from each station element can also be ana- 
lyzed with real-time logic and/or copied to LOFAR' s Tran- 
sient Buffer Boards (TBBs), circular RAM buffers holding 
up to 1.3 seconds of data from each station element. Data 
from the TBBs can be accessed and returned for processing 
without interfering with other observation modes. 

The LOFAR Cosmic Rays Key Science Project is using 
this parallel capability to search for atmospheric cosmic 
ray events. The LBAs individually provide all-sky cover- 
age, and their sensitivity peaks in the 30-80 MHz range 
where the maximum radio signal from extensive air show- 
ers (FAS) has been observed [3]. The LOFAR "VHECR" 
(very-high-energy cosmic ray) mode utilizes the LBAs for 



a real-time search for FAS signatures, and its implementa- 
tion and operation is the subject of this contribution. The 
"HECR" (high-energy cosmic ray) mode, which trades sky 
coverage for increased sensitivity by using station beams to 
search for FAS, will be developed in the near future, while 
the "UHFP" (ultra-high-energy particles) mode, by which 
beams from all stations target the Moon to detect particle 
interactions in the Lunar regolith, is the subject of another 
contribution [7]. 

2.1 Triggering and Data Acquisition 

In VHFCR mode, LOFAR can trigger on cosmic rays in 
two ways: directly on the radio air shower signal, or via an 
external trigger (e.g. from a particle detector array). Radio 
triggering is performed in two stages, and the time-series 
signal around a pulse is stored for all channels once the 
highest-level trigger is satisfied. 

Our first-level trigger is based upon that developed for the 
LOPES experiment [9]. Pulses are detected in real-time, by 
processing the digitized radio signal of each polarization 
channel in an FPGA (field-programmable gate array). The 
FPGA performs a threshold search on the absolute values 
of the voltage ADC samples using the criterion 

\xi\ > k^ii, (1) 

where /i^ is a running average over the previous 4096 ab- 
solute values, and k is an integer threshold factor, typically 
set between 6 to 8. If we assume the background signal 
to be Gaussian-distributed with standard deviation a, the 
mean absolute value is related to a as 




and thus the threshold is typically a factor of 5 to 6 
above the time-varying noise level. When this threshold 
is crossed, a message is passed to the station's control PC, 
where the second-level trigger operates. The message con- 
tains a number of parameters characterizing the pulse, in- 
cluding trigger time, pulse height, pulse width, and noise 
levels before and after the pulse. 

The second-level trigger searches for coincidences between 
single-channel triggers. Typically, we require 16 to 32 
out of 96 channels to trigger within a 1 /is time window. 
For these events, we perform a real-time direction estimate 
based on the first-level trigger times and the antenna po- 
sitions. Although the uncertainty in time of the trigger 
messages is about 2 to 3 samples (i.e. 10-15 ns), the es- 
timated source direction is accurate enough to distinguish 
between sources from the sky and interference sources on 
the horizon. The coincidence event rate is dominated by 
a few fixed radio-frequency interference (RFI) sources — 
most commonly electric fences around nearby farm fields 
— that emit pulsed radiation. Therefore, the second-level 
trigger also uses this preliminary reconstruction to exclude 
triggers coming from within 30° of the horizon. 

L A 160 MHz mode is also possible, but is used less frequently. 
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When a pulse passes all second-level trigger criteria, the 
buffered time-series data for 0.5 to 1 ms around the pulse 
are stored for offline analysis. The relatively large size of 
this data window allows offline spectral cleaning of narrow- 
band radio transmitters from the signal with a frequency 
resolution of 1 to 2 kHz. 

In addition to the radio-triggering mode, LOFAR can also 
be triggered by an external particle detector array. LORA 
(the LOFAR Radboud Air Shower Array) is an array of 20 
plastic scintillator detectors arranged in 5 stations within 
the LOFAR superterp [8]. Upon detection of an air shower 
by LORA, a trigger can be sent to LOFAR to store the radio 
signal from all dipoles for further analysis. This can be 
used to provide an unambiguous cross-check that a radio 
signal is from an EAS, and can also be used to reduce the 
detection energy threshold, as described in Sec. 2.3. 

2.2 Data Analysis 

We process the second-level trigger data using an auto- 
mated pipeline structure, the first objective being to dis- 
criminate cosmic ray (CR) events from other pulses. The 
following steps are performed: 

1. data vaUdation, e.g. elimination of pulse trains or 
other noise; 

2. matching of the event data with coincident triggers 
from particle detectors and/or other LOFAR stations; 

3. spectral cleaning, i.e. removal of narrow-band radio 
transmitters from the spectrum; and 

4. reconstruction of the source position (direction and 
distance). 

To locate the source, we calculate the "beam-formed" time- 
series for a set of directions and distances. Assuming a 
point source at a given position, we correct for the relative 
time delays at each anterma, using a spherical-wavefront 
approximation, and then sum the delay-corrected signals 
from all antennas. As the true source position can be iden- 
tified as the one yielding the maximum beam-formed pulse 
power, the best-fit direction can be efficiently found using 
a downhill-simplex algorithm, using the second-level trig- 
ger's direction estimate as a starting point for the search. 
The results of the pipeline can then be used to determine 
key air shower event parameters such as energy and posi- 
tion of shower maximum. 

2.3 Energy Threshold and Estimated Event Rate 

The detection energy threshold of the VHECR mode is de- 
termined at minimum by the first-level trigger on an in- 
dividual LBA dipole. An event with signal exactly equal 
to the threshold will trigger on average half of the dipoles 
at station level (thus passing the multi-channel trigger re- 
quirement), and will not be rejected as RFl unless it is 
a highly-incUned event with zenith angle larger than 60°. 
Thus, the second-level trigger restricts the sky to tt sr for 
real events. 

The noise background against which cosmic ray signals 
must be detected is dominated by the Galactic background 



emission, with effective temperature at wavelength A given 
by Tsky = 60 ± 20(A/m)2-55 K [2]. The sky tempera- 
ture is in turn defined using the measured flux via 
the Rayleigh-Jeans Law of = 2fcBTskyA~^, where 
ks is Boltzmann's constant. The (sparse) outer LBAs 
have an effective area Adi that increases with wavelength: 
Aeff « A^/S. As there is strong RFI below 30 MHz, we 
use a filter which Umits our bandpass to approximately SO- 
SO MHz; after filtering, RFI typically contributes an addi- 
tional 25% to our system noise. Assuming constant effec- 
tive area over the observed solid angle Osky and account- 
ing for RFI, the noise power P„ received in a single polar- 
ization channel can be calculated by integrating over fre- 
quency v. 

Pn = 1.25 Q,ky j ^^^^ A,s{v) m du (3) 

where f{u) is the power response of the filter. The radio 
signal Vi^as of an extensive air shower is coherent and must 
be calculated in the voltage domain from its intrinsic field 
strength 5(2/) and angle 6p of the receiver to the polariza- 
tion vector: 



Feas = COS Op J S{v) f (ly) A^e (i^) dv. (4) 

This must sufficiently exceed the RMS noise voltage 
Kms — PnZo in order to result in a trigger. Simula- 
tions of the expected signal S{p) indicate an increased 
strength at lower frequencies [4]; however, we use a em- 
pirical frequency-independent parameterization based on 
LOPES measurements [10]. 

Using a trigger threshold of |ycas| > 5|\4ms|> we find a 
minimum energy threshold of ~ 5 x 10^^ eV. However, 
the strength of the geomagnetic emission mechanism is 
strongly dependent on the angle between the shower axis 
and the Earth's magnetic field, and a typical shower geom- 
etry results in a threshold closer to 10^^ eV. Taking the sta- 
tion effective area to be equal to that of a circle with radius 
given by the characteristic signal fall-off length of ~ 200 m 
[10], the total expected event rate for all 40 core and remote 
stations will be of the order of 1/hr. The coincidence rate 
of the superterp stations with the LORA particle detector is 
estimated to be of order 1/day. 

For LORA-triggered events (where the arrival direction is 
given by the particle detectors), the radio threshold is given 
instead by the beam-formed data from all 48 antennas. 
Thus the threshold for LORA events will be lower than 
radio-only mode by a factor of between and v^, de- 
pending on the relative strengths of the signal in each of the 
polarization channels (2 per antenna). Therefore we expect 
a threshold of ~ 10^^ eV over the central region covered 
by LORA, and a corresponding event rate of order 2/hr. 

2.4 First Cosmic Ray Events and Outlook 

The first confirmed air shower radio pulses were observed 
wifli LOFAR in June 2011, using LORA as a trigger. The 
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Figure 2: A cosmic ray air shower observed with five LO- 
FAR stations (groups of circles) and LORA (pentagons). 
The circle size represents the power of the radio signal, and 
the color scale indicates the arrival time. The reconstructed 
shower core and arrival direction are shown by the cross 
and dotted line. 
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Figure 3: Time-series air shower radio pulses from LO- 
FAR antennas (superimposed), after correction for geomet- 
ric arrival time delays. 

cosmic ray nature of the pulses is confirmed by compar- 
ing independent directional reconstructions by LOFAR and 
LORA, which typically agree to within a few degrees. 

A sample air shower event is shown in Figs. 2, 3, and 4. 
Fig. 2 shows the time delay and signal strength across the 
core stations. Radio pulses from individual dipoles ex- 
hibit coherence after correction for geometric delays (see 
Fig. 3). A preliminary lateral distribution function of the 
same event, showing the radio pulse power as a function of 
distance to the shower axis for the two LBA polarizations, 
is shown in Fig. 4. The error bars are preliminary and are 
based on average LORA reconstruction errors. The rela- 
tive normalization between the two polarizations does not 
yet take into account differences in directivity, and a full 
calibration of the radio signal into physical units (electric 
field) is still in progress. 




100 150 200 

Distance to Shower Axis [m] 



Figure 4: Radio pulse power as a function of distance to the 
shower axis, for a single air shower event seen in five LO- 
FAR stations Both LBA polarizations are shown (aligned 
northeast-southwest and northwest-southeast). 

These data are, to our knowledge, the most densely instru- 
mented measurements of radio air shower emission to date. 
Analysis of the lateral distributions is underway and will 
allow detailed studies of the polarization and falloff char- 
acteristics of the radio signal, including effects from sub- 
dominant emission mechanisms. Studies of the wavefront 
shape are also planned. 

In addition to the LORA-triggered mode, optimization of 
the radio self-trigger mode is continuing. The improve- 
ments underway will allow a reduction of deadtime as well 
as the veto of RFI sources from above the horizon, such 
as those from airplanes. The self-trigger mode will signif- 
icantly increase the detection area available for cosmic ray 
observation. 
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